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Abstract 

Symmetry energy terms from macroscopic mass formulae are investigated as generalized polar- 
izabilities of nuclear matter. Besides the neutron-proton (n-p) symmetry energy the spin dependent 
symmetry energies and a scalar one are also defined. They depend on the nuclear densities (p), 
neutron-proton asymmetry (6) , temperature (T) and exchanged energy and momentum (q) . Based 
on a standard expression for the generalized polarizabilities, a differential equation is proposed to 
constrain the dependence of the symmetry energy on the n-p asymmetry and on the density. Some 
solutions are discussed. The q-dependence (zero frequence) of the symmetry energy coefficients 
with Skyrme-type forces is investigated in the four channels of the particle-hole interaction. Spin 
dependent symmetry energies are also investigated indicating much stronger differences in behavior 
with q for each Skyrme force than the results for the neutron-proton one. 
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1 Introduction 

The symmetry energy terms and their dependences on the density are of relevance for the nuclear 
structure and in many nuclear processes including the structure and dynamics of proto-neutron and 
neutron stars. The neutron-proton symmetry energy is the best known in spite of the different val- 
ues at the saturation density in the literature (in the range of 25MeV up to 36MeV). It is basically 
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represented by a squared power of the neutron-proton (number or density) asymmetry in usual macro- 
scopic/microscopic mass formula the parabolic approximation. With a symmetry energy coefficient 
(s.e.c), a-r, the binding energy is, in the simplified versions, usually written as: 

E/A = Ho{A,Z)/A + ariN - Zf/A'^, (1) 

where the energy density Hq does not depend on the asymmetry, Z, N and A are the proton, neutron 
and mass numbers respectively. Neutrons and protons occupying the same total volume yield a term 
proportional to the squared asymmetry density, ar{pN—pz)'^/p^ which appears in the nucleonic matter 
equation of state. The neutron and proton densities may not be exactly equal to each other in nuclei 
Different polynomial terms of the asymmetry in this expression (proportional to {N — Z)" for 
n / 2) are usually expected to be less relevant [Q, ^, ^, |5|, |^, 0, ^. However it is not well known 
whether and how this parabolic approximation is to be modified for very asymmetric systems, such as 
nuclei far from the stability line or for (asymmetric) nuclear matter above and below the saturation 
density In large stable nuclei such as "^^^Pb the n-p asymmetry {{N — Z')^ jA? ~ 1/9) is not so 
large as it would be in neutron matter. The n-p symmetry energy coefficient (s.e.c.) is given by 
the static polar izability of the system [11^] (the inverse of the "isovector screening function") which 



can also depend on the asymmetry of the medium |11|. This may lead to slightly different forms for 
the symmetry energy for very asymmetric n-p systems. 

Other symmetry energy coefficients may also be defined in nuclear matter, for instance, the spin 
Afj and spin-isovector A^t ones. Extending the n-p symmetry energy, the other symmetry energy 
coefficients can be defined in macroscopic mass formulae as: 

_ -ffo(^, Z) ^ ^ {pN — PzY _|_ ^ {Sup — Sdownf _|_ ^ {Pup ~ Pdown + Pdown ~ Pup?' f^c^^ 
A P ^ {Sup + SdownY P^ 

where the density (and eventually number) of neutrons and protons is denoted by Pn^Pz^ of nucleons 
with spin up (down) by Sup{Sdown) and Pupdown neutron/proton densities with spin (up, down). 
The spin channel may lead to the appearance of polarized nucleonic matter which has been investigated 
within different approaches with controversial results |12, 13, 14, 15, 16, |ll], [l^, 19, 2C, 21|. The 



spin channel is also relevant for the study of the neutrino interaction with matter because it couples 



to the axial vector current together with the scalar channel in dense stars 23, 24 1. The spin- 
isospin channel has been associated to pion condensation p5| , ^ and also to anti-ferromagnetic states 
[20|. A nuclear dipolar incompressibility was also defined in being related to the nuclear matter 



incompressibility as discussed below, and which varies accordingly with the n-p asymmetry being 



eventually relevant for the isoscalar dipole resonances |27]. These coefficients and their corresponding 
dependences on the asymmetry of neutron-proton densities have been investigated in several other 



works, as for example in |28, 29, 21]. A different way of obtaining the symmetry energy has been 
proposed by means of the linear response method for the dynamical polarizabilities. The static limit 
of these generalized polarizabilities are proportional to the inverse of the symmetry energy coefficients 
in symmetric matter ||l^. Therefore it becomes reasonable to consider the polarizabilities as a suitable 
and sound framework to determine the behavior of the symmetry energy with the parameters of the 
nuclear equation of state. Developments with relativistic models also yield strong effects with the 
isovector mesons, see for example references |30, ^ among others. The density dependence of 
the neutron-proton symmetry energy and the isospin dependence of the nuclear equation of state are 
being extensively investigated for several reasons and several experimental tests are being done and 
prepared nowadays mainly in intermediary and high energy heavy ion collisions |33, 34, 35, 36, 38 



p9| , 41, 42 1 . The investigation of the possible effects with their particularities and the consequences 
for the observables is extense involving new experimental facilities such as RIA and GSI |36, ^ 
besides many other works |^^. Any definitive realistic investigation at really high densities (several 
times the saturation density) should take into account baryonic structure with the internal (quark and 
gluon) 

In the present work some aspects of the symmetry energy terms are investigated as provenient from 
the generalized polarizabilities of nuclear matter for different ranges of the density, n-p asymmetry 
and momentum exchange, in the zero energy limit of the dynamical polarizabilities, within general ar- 
guments and with Skyrme forces at finite temperature. The case in which there is also non zero energy 
exchange corresponds to the analysis of the dynamical response function. The density dependence of 
the equation of state is not well known and it is reasonable to ask whether and how the n-p symmetry 
energy (and more generally other symmetry energies in the other channels of the nuclear interaction) 
depends on isospin at different densities / very high n-p asymmetries. The parabolic approximation, 
usually appropriated for a restricted range of densities (very) close to the saturation po and small 
asymmetries, may be modified for lower and/or higher densities. 

The parameters of the forces which are used (SLyb and SKM) were fitted (i) from results of 
asymmetric nuclear matter and neutron matter properties obtained from microscopic calculations 



[46 1 and (ii) from properties of giant collective modes in ^"°Pb [47|. Other forces will be investigated 
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elsewhere. Skyrme forces can be obtained from a reduction of the nuclear density matrix [48| and their 
basic structure is also present in non relativistic reductions of relativistic models for nuclear systems by 
passing to relativistic point coupling models or not |p9|, pO| such that the necessary density dependence 



of each of the terms are expected to be stronger than considered in the earlier parametrizations |51, |52| . 

This work is, in part, an extension of previous works and it is organized as follows. In the 
next section general aspects for the investigation of symmetry energy within the approach of the 
general polarizabities are discussed including the stability of nucleonic matter with respect to external 
perturbations. In section 3 an analysis of simultaneous dependence of the polarizabilities on the 
neutron-proton asymmetry and on the density is proposed with a differential equation that constrains 
these two behaviors of the symmetry energies. In sections 4 and 5 the g-dependence (exchanged 
momentum between the components of nuclear matter, eg. neutrons and protons) of static generalized 
polarizabilities at finite temperatures with Skyrme forces is investigated in the limit of symmetric 
nuclear matter. In the last section results are summarized. 

2 Symmetry energy and nuclear matter polarizabities 

Basically, in this section, arguments from previous works are reproduced. Consider that with the 
inclusion of an external source of amplitude e, which separates nucleon densities with quantum numbers 
(s, t) (where (0, 1) is for spin up-spin down and (0, 1) for neutron-proton), the energy density of nuclear 
matter can be written as: 

H = Ho + A JlMl^fMll + (3) 
P 

where Hq does not depend on the density asymmetry {p(^s,t)i ~ P(s,t)2)^' -^s,* corresponding 
symmetry coefficient (^i,o the spin one, Ao^i the neutron-proton one) and the total density fluctuation 
is (3 = 5p(^s,t)i ~ ^P{s,t)2i these two cases. For the spin-isospin external perturbation {s, t = 1,1) the 
simultaneous fluctuations of the spin (up/down) and neutron/proton densities are to be considered 
just like it is shown in expression (|2|). In the case of (s = 0,t = 0), the scalar channel, there is a 



change in the total nuclear density and Aq^ is associated to a dipolar incompressibility [11 1. For equal 
volumes the densities become the nucleon-numbers. 

In the ground state, the variation of the energy with respect to the density fluctuation of a channel 
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(s,t), 6p = /?, yields the condition of minimum: 

e' + 2^(p„ + M = e + 2^6p= 0, (4) 

where pm = Pq — Pq ^ is for an n-p asymmetric matter (or correspondingly pf^ = p^ — pg°"'" ^ for 
spin polarized matter) and the total "inducing perturbation" for n-p asymmetric systems is denoted 

by 

e = e' + 2As,t Pm- (5) 

The ground state can be considered to have a (polarized) spin up-down asymmetric density given by 
7^ simultaneously to (or instead of) the n-p asymmetry (also denoted by pm above), pm will be 
considered in most part of this paper . The nuclear matter polarizabity in the channel (s, t) can be 
written for e' or for the total (inducing) perturbation, e, respectively as: 

nfEE^ = — P—, ii) w^'^^ = — {ti). (6) 

The stability condition for these expressions are different. This will be discussed below. 

The main development will be focused for the neutron-proton symmetry energy (s, i = 0, 1) al- 
though it is analogous for the other channels. The neutron proton asymmetry used in the present 
work is defined by the neutron and proton densities Pn, Pp as: 

6=^-1. (7) 

An asymmetry coefficient which is probably more familiar to the reader is given by: 

a = ^^P--P\ (8) 
P 

They are related by: b = 2a/(l — a). The coefficient b varies from 6 = 0, in symmetric nuclear matter, 
up to 6 ^ oo, in neutron matter. For the sake of generality the coefficient As^t is considered to be a 
function of the density fluctuation (3. The fluctuation (3 is considered to depend on the n-p asymmetry 
b. These parameters may be related to each other and therefore it will be written that Ag^t = Ag^tiP)- 
Relations between b and (3 have been investigated by means of prescriptions. Among those, one which 
leads to reasonable results is: 

^ = ^Pn{P^^, (9) 



,1 + 6. 

Where 5pn is the neutron density fluctuation. In the n-p symmetric limit (3 = 26pn and in another limit, 
in neutron matter, (3 = 5pn- This ansatz (expression ^) is based on the assumption that the density 



fluctuations are proportional to the respective density of neutrons and protons, i.e., 5pn/ = Pnl P-, 
being p the total density. 

The resulting expression for the symmetry energy coefficient Ag^t for the prescription above is given 



by inl: 



for a general s, t (spin,isospin) channel of the effective interaction. In this expression Agym = a-r — 
30MeV is the s.e.c. of symmetric nuclear matter (6 = 0). The (generalized) n-p symmetry energy 
term can be rewriten as: 

A ' = Agyjji{l + a), 

which corresponds to a third order term in the binding energy, being smaller than the quadratic term 
because q < 1. In not very n-p asymmetric systems, those with n-p asymmetry close to the stability 
line, << 1. For b = 2 (a = 0.5, neutron density three times larger than the proton density) 
it follows A = 1.5Asym- In the limit of neutron matter A{b oo) = 2Asym- For proton excess 
6 < 0. Prescription (|9|) is therefore model-dependent and different choices yield other forms for the 
the (asymmetric) static generalized "screening functions". The dynamical response functions are less 
sensitive to this prescription. 

So far it has been assumed that the stable density p is independent of b (or a). Below it is envisaged 
a development to guide the simultaneous variations of these two variables. Several investigations of 
the role of the symmetry energy on observables in Radioactive Ions are being prepared for RIA and 
GSI. For this the density dependence of the symmetry energy is extremely relevant. However for 
mass formulae of very asymmetric nuclei and for the equation of state at densities different from po 
the isospin dependence of the symmetry energy may be different from the usual one given by the 
parabolic approximation, expression (|l|). Furthermore, more elaborated pictures in relativistic mean 
field calculations, which considers isovector mesons, 6, yield a qualitative increase of the relevance 
of the neutron-proton asymmetry, with a larger difference of neutron and proton effective masses 



[54, 55, 32]. Experimental bounds on the neutron and proton effective masses |56] may shed light 
on this. The spin and spin-isospin symmetry energies can be investigated analogously. For example, 
the behavior of ^i^o of the spin channel (which has already been written as in the static limit in 
the framework of the Landau's Fermi liquid theory), at variable densities was investigated in different 



works [13, |13|. 
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2.1 Stability conditions 

In the usual case in which the density p is not dependent neither on b nor on (3 there are two ways 
of writing a solution for the polarizability Yis,t from expression (^). They correspond to the different 
definitions of the external source shown before, respectively e' (i) and e = e' + pm They allow for 
defining polarizabilities given respectively by: 

where C^*'*^ = — constant, with the usual value of the symmetry energy coefficient. In 

the n-p channel: (-4g^™)" ~ 30 MeV (symmetric limit). These two polarizabilities (11) are equal in 
the limit of symmetric nucleonic matter pm = 0. This derivation applies for any of the channels (s, t). 

Consider that the binding energy is to be minimized with respect to the density fluctuation (3. From 
this an equilibrium condition for nuclear matter is obtained with 5"^ {E / A) / 5 13"^ > 0, being different 
from other ones and complementary to them |^]. To be a stable minimum of the binding energy the 
coefficients, of both definitions of the polarizabilities, satisfy respectively: 
6'^E/A 2Cst , , S'^E/A 

^ = -(n5y>0. (')■ ^ = 2A.>0. (n). ,12, 

where 11" ^ and Ag^t given in expressions ([Tl| (i)) and (|ll| (ii)) respectively. The constant C*^*'*^ may 
be negative (stable symmetric nuclear matter pT| ) or positive, As^t and Ilg^t also may be negative or 
positive. While the second expression yields the more expected result, i.e. the stability is directly 
shown by the signal of Ag^t in each channel {s,t), the first expression has a more involved behavior 
due to the complicated form of expression ( pT| (i)). From the polarizability (|Tl|-(i)) two conditions for 
real A and stable system follow: 



(i) (1) Al, < {Al,)sym, (i) (2) {Al,)sym {2 - ± 2. /l - f'' ] < 0. (13) 

\ {•^s,t)syin y {■^s,t)syrn ) 

From condition ( [T3| (i)-(l)) the neutron-proton asymmetry can only lower the value of the generalized 
coefficient Ag^t to keep the system stable with the use of polarizability ([l^ (i)). It is worth emphasizing 



that the two conditions (12) with the respective definitions for Il^^t should not be mixed. If the 
polarizability from expression ([Tll-(ii)) is considered the condition ([T2|-(ii)) is to be applied, otherwise 



inconsistent results arise. Expression (11- (ii)) is the usual form. However if one considers solution 
(|Tl|-(i)) the condition ([l^(i)) is to be applied, otherwise inconsistent results arise. In particular in 
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the case of the polarizabity given by expression (|Tl|-(i)) there are several possibihties for the stabihty 
of the symmetric and the corresponding asymmetric matter depending on As,t > (or As,t < 0) and 
■^sym > (o^ •^sym < 0) each of the channel (s, t). The microscopic in medium nucleon interactions, 
in an exact calculation, would give the correct one. Expression ([lO| ) for Ao,i{b) was found with the 



solution (11 (ii)) for Hg^t- 

The stability conditions of a Fermi liquid in the leading order, in each channel of the interaction, 
correspond to a particular case of the above expression ( [T2| - (ii)). They are given by the denominator 
of a particular limit of the response function Ils,t which can be written as: 

as,t = No{l + Jo''*) > 0, (14) 

where Jg'* stands for any of Fq, Fq, Gq, G'q respectively for the scalar (s = 0, t = 0), isovector (s = 0, t = 
1 with Ur), spin (s = l,t = with a^) and spin-isovector (s = l,t = 1 with a^-T-) channels [E^, 59|. 



These expressions contain the leading terms of the more general calculation. Within a non relativistic 
formalism with Skyrme type interactions they can be written in terms of Landau parameters ||l^, |60| . 
Other considerations can be associated in different formalisms pT| , p^ , |l3|, |6^. A complementary and 
more general discussion for particular models will be done in a forthcoming work. 



3 Simultaneous dependence on isospin and density 

Next it will be assumed that there is an implicit and a priori unknwon dependence of the saturation 
density on the n-p asymmetry without any supposition about the microscopic origin for this, po = 
Po{b). From the general and usual expression for the polarizability (|6|) (or ( p!l| (ii))) a differential 
equation for the simultaneous isospin and density dependence of the symmetry energy (coefficients) 
As^t will be derived. Although expression (^) was also derived without considering a dependence of p 
on b it will be considered that this simple form is more general. The derivative of the polarizability 
Ils,t, expressions (^, with respect to b is given by: 

dp _ 1 dAsA dp 1 dAs,t] . . 

db ^\[p As,t dp J db As,t db j- ^ ' 

The variation 5 (5 /5b is given by expression (^), the prescription for the relation between the fluctua- 
tions. This equation has other three derivatives a priori unknown which have to be consistent with 
the equation of state: the derivatives dA/dp, dA/db and dp/db. This expression is therefore to be 
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equated to that of prescription (0) or, more generally,: 



f = -mb). (16) 



The resulting equation is: 



1 1 OAs A dp 1 dAs,t 



-fib)- (17) 



[\p As,t dp J db As,t Ob j 

This is one of the most relevant results of this paper. This equation constrains the simultaneous 

dependence of the symmetry energy on the density and on the nucleon density asymmetry (through 

the generalized coefficient As,t, in the channel s,t, which is not a constant anymore) 0. 

The following cases correspond to the derivation of section 2 (expressions (pA])): 

dp dA A 

— = 0, and / or ^ = — • 
ob op p 

These correspond to the use of prescription given by expression (^) which yields the function: 

For this prescription, which yields expression ([lO|) for A{b), it has been assumed that p was independent 
of b, therefore, in that case: 

= (19) 
Pdb Adb' ^ ' 

In this case the behavior of Aip) can be the one typical of relativistic models with the increase of (any 
of the) symmetry energy coefficient with the increase of the nuclear density, i.e., As,t oc p and in part 



of microscopic approaches, for which its value usually tends to a constant, |2^, |6J, ^ |6g, p7 |. 
However this is not the most general and interesting case because condition ([l^ holds when p is 
independent on b. 

A slightly more general parametrization can be investigated. For example, Heiselberg and Hjorth- 
Jensen []68| , |69t| have used the following expression for the density dependence of the symmetry energy 
which nearly summarizes results obtained from relativist models: 

Esym = EsymiPo) y~ J ' (^'-*) 



^The parameter b however may be replaced by the equivalent one for the spin up-spin down asymmetry in polarized 
nuclear matter. The same form is obtained for spin-polarized nuclear matter, by interchanging the neutron-proton 
variables (from s,t = 0, 1) to spin up- spin down ones (with s,t = 1,0). An equivalent prescription has to be provided 
for the spin density fluctuations. 
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where 7 is a constant and po the saturation density. A variational calculation favor values of the order 



of 7 ~ 0.6 whereas an analysis of heavy ion collisions experiments at low energies 7 ~ 2 [68, 69|. Prom 
the differential equation ( [TtI) it will be considered parametrizations given by: 



dA_ A 

dp ~ 

A = A. 



7 

P 



9P , (21) 



as + a^b 

where Ui (1=1,2,3,4) are constants. When the asymmetry coefficient in expression ( p^ reaches the 
value 

the symmetry energy coefficient ^0,1 changes the sign making the system unstable according to the 
stability condition (|l^ (ii))- The resulting equation, from the general equation ([ITD, for the density as 
a function of h with the above parametrizations is: 

For positive as and the general solution is given by: 

p{h) = {h + 2f~> (1 + hY~~^ (as + a4 bf ' e^-"'^' B. (23) 

In the limit of symmetric nucleonic matter, for b the neutron-proton density asymmetry, 6 = 0, the 
constant B can be fixed in terms of pQ. In neutron matter p ^ or p ^ cx3. For this to be finite, 
p{b 00) — > 0, one must have a2/a4 < and aia4 — a2as < 0. A particular solution appears for 
03/04 = —b which yields p = 0. When this occurs ai and a2 from expression ( |23| ) have different signs. 

For the usual form for the symmetry energy term in which ^ = a,- is independent of 6, the resulting 
density as a function of the asymmetry b from the differential equation (^) is given by: 

Pib)=cJp^)\ (24) 



,1 + 6. 

where Cq is fixed by a boundary condition, for example p{b = 0), with a fixed value of 7. From this 
limit: Co = po2^^'^ whereas in neutron matter p{b 00) = Co- The ratio of the density in these two 
limits is given : 

,(6^oo) = ^(^^. (25) 

2-r 



However, the parameter 7 from the parametrization ( pO| ) may be (assumed to) depend on the neutron- 
proton asymmetry coefficient b (or equivalently a). In this sense, a modification in the usual symmetry 
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energy dependence on the n-p asymmetry can be expected to be equivalent to different values for the 
parameter 7, at different densities, in different experimental situations. 

A different form for the above equation ([T7| ) can be written by considering that: ^ = g{b,p), and 
^ = h{p, b) ^ ^. The following differential equation appear for A{b, p) with these functions: 

(20) 

Considering the particular prescription (|T^ ) it is obtained the following expression: 

(^-^- -^h{P,b)^ 9{b,p) = -2f{b). (27) 

These expressions can be considered for any channel {s, t). They generate one differential equation 
for each channel of the nuclear effective interaction with Ag^t, and therefore the final p dependence on b 
is to be the same for each of these equations, for b representing the same asymmetry (neutron-proton, 
spin- up-spin-down). For this, the choices for As,t and /3(6) should be associated, otherwise there will 
appear different p{b). 



4 Generalized "Screening functions" with Skyrme forces 



In this section the analysis done previously [10, 53, 11, 21] is extended with the static limit of the 
expression for the dynamical polarizability of a non relativistic hot asymmetric nuclear matter with 
Skyrme effective interactions, limuj-,o^s,ti^, q)- These polarizabilities were obtained by the calculation 
of the response function of hot asymmetric nuclear matter in terms of three densities: neutron and 
proton densities {pi), momentum density (rj) and kinetic energy (jj) densities from the time dependent 



Hartree Fock approximation with Skyrme forces [10, 11 [. These densities appear in reductions from 
relativistic models in which the scalar density is written in terms of them [4^, 5C]. The time dependent 
approach introduces CP violating terms proportional to j which are larger in asymmetric nuclear 
matter. Four asymmetry coefficients are defined, a, b, c and d for the effective masses and densities, 
and they are given by: 

^ _ Pon , 1 + ft , 1 

Pop 



1, 



1, 



d 



(28) 



Pop 2 + 6' 1 + (1 + 5)3 

where m* are the neutron and proton effective masses. Small approximations were done: (i) to equate 
the asymmetry coefficient defined for the momentum density to the density asymmetry coefficient (ii) 



11 



to choose a particular prescription for the fluctuations of the asymmetry density - expression @). The 
second approximation is in fact a choice with dynamical content and it deserves more attention. 

At the Hartree Fock level, the symmetry energy coefficient, ar = O'^^Eq/ A)da^ , can be written 
from the expansion: 



E, d{EQ/A) 
— = Ho + a 



A ^ da 



d\Eo/A) 
2 da^ 



+ ..., (29) 

a=0 



a=0 

where higher order terms are not written. There may appear (small) higher order terms. By calculating 
the general polarizability, within the linear response approach, a whole class of ring diagrams contribute 
beyond the Hartree- Fock |7^. Therefore more complete symmetry energy terms can be obtained. 

For the calculation of the response function an external source is introduced in the Hartree Fock 
time-dependent equation which induces small amplitude density fluctuations. The general form of the 
source is the plane wave one, it is given by: 

V,xt = -eds,tVe-'(-'-^-^\ (30) 

with an amplitude e (usually a small parameter), an associated dipole moment D (equal to the unit 
from here on) and the operator O acts on the nucleon states. In particular, for the isovector interaction 
the third component of the isospin (Pauli) matrices is considered yielding neutron-proton density 
fluctuations. With the above external source, the Hartree Fock equation for nuclear matter is written 
as: 

dtPi = -i[Wi + vr\pi], (31) 

where Wi is the Hartree- Fock energy of protons or neutrons. The induced density fluctuations Sp are 
to have the same spatial and temporal plane waves behavior of the external source. 

The resulting expression is more appropriatedly written in terms of generalized Lindhard functions 
whose real parts, at zero temperature F2i, were defined as fl^, |53|: 



M,^(..,q) = -^^ej d fe ,^,^_,.(i,)^,.(i,^q) (k.(k + q)) =^J d^me F,.. 

(32) 

In these expressions fi{k) are the fermion occupation numbers for neutrons (i = n) and protons {i = p) 
which will be considered only for the zero temperature limit (when dfi{k) —5{k — kp)), g is the 
degeneracy factor for spin and isospin, M* is the effective mass in symmetric nuclear matter. In the 
limit of zero energy exchange {iv — > 0) the Lindhard functions yield the (q-dependent) proton and 
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neutron densities, momentum and kinetic energy densities are given by: 

7M* 



Pi = ^ l dfiik) ^eF^iiv ^ 0), (33) 
Mq = ^ j dfiik) ^eF^iu ^ 0). 

In the symmetric nuclear matter the momentum-dependent polarizabihty (P) in the channel s, t 
for the Skyrme effective force parametrization is written as: 

AsAq) = ^ {1 + '^Wn" + ey/'^MV? + iV,''')^iM*)^ (dp^^ - AM^Ng) } , (34) 

Where Vb(g^) and Vi are functions of the Skyrme forces parameters for each of the (s,t) channel 
shown in the Appendix. The nuclear matter incompressibility modulus is related to Aofiiq^ = 0) in 
the Appendix. The q-dependent densities N'^, p^, are the total densities from expressions (|33|). The 
term proportional to can be re-written in a homogeneous nuclear matter at zero temperature, as 
pr—j^ which is to be zero in the Galilean invariant (homogeneous and static) limit |5C]. This invariance 
is broken in these cases and it is amplified in asymmetric nuclear matter, q in the neutron-proton 
channel is the exchanged momentum between the neutron and proton components, and similarly, 
in the spin channel, the corresponding exchanged momentum for spin up and down nucleons. The 
stability condition for this expression is given by (^ - (ii)). The finite temperature calculation of 
the densities lead to finite temperature symmetry energy coefficients. In the zero frequence limit the 
imaginary part of the response function disappears. 

In the limit of low momenta, q « 2kF, the w = limit of the Lindhard functions are simplified, 
as shown in the Appendix. The polarizabilities of symmetric nuclear can be approximatedly written 
in the following form: 

As,t = AsAT, P) + A^S P)9 + <i P)l'^ (35) 
where Ag^T^p) is the usual symmetry energy coefficient in the channel (s,t) |l|] and A^*|(T, p) 
are functions of the Skyrme force parameters (combined in the functions V^'^ and T^{*'*), p and T. 

4.1 Results for Skyrme interactions in the 4 channels 

In this section the generalized polarizabilities are investigated numerically for the Skyrme forces SKM 
and SLy(b) for the four channels of the particle-hole interaction as functions of the exchanged momen- 
tum at the normal density po- For this, the chemical potential was adjusted to mantain a constant 
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stable nuclear density p{T) = po- As a consequence the results are not very strongly dependent on T. 
In Figure 1 the neutron-proton polarizability is shown as a function of the (exchanged) momentum for 
temperatures T = 0, 4 and 7 MeV for Skyrme force SLyb and at T = MeV for the force SKM. Both 
forces produce widely accepted values for the symmetry energy coefficient, Ao^i{q = 0, T = 0) ~ 32 
MeV. There is a general behavior (for both forces at any of the temperatures) of decreasing Ao^i with 
increasing exchanged momentum up to g ~ 500 MeV. This corresponds to nearly twice the nucleon 
momentum at the Fermi surface. However ^o,i does not reach negative values. The behavior of de- 
creasing values of A^'^ for increasing momenta is in agreement with other analysis for the momentum 
dependence of the symmetry energy [^]. The symmetry energy coefficient has, according to expres- 
sions (3^,35), a linear/quadratic behavior for very small exchanged momenta q. This is followed by an 
abrupt change of behavior at g ~ 500 MeV. For higher q the generalized s.e.c. (^o,i) increases nearly 
linearly at different temperatures. For the s.e.c. Aq^i to become negative the Skyrme parameters Iq 
and should result in larger values of Vq than those of SLyb (this variable is still smaller for the 
force SKM) and/or different values for ti,t2- 

In Figure 2 the spin-isospin generalized polarizability dependence on exchanged momentum be- 
tween neutrons and protons with spin up and down is investigated for the same cases of Figure 1. 
There is again a quite defined change of behavior at g ~ 500 MeV. The generalized spin-isospin s.e.c. 
remains nearly constant with increasing g up to g ~ 500MeV. The force SKM yields smoother varia- 
tions than SLyb like in the n-p channel. Above the saturation density, ^i^i decreases for most forces 



eventually reaching a negative value |21| . Finite temperature effects are larger for higher values of q. 

In Figure 3 the spin generalized polarizability, Ai^, is shown for the same cases of the previous 
figures. The turning point present in the isospin-dependent channels, investigated in figures 1 and 
2, is the same {q ~ 500 MeV). This is due to the form of the Lindhard functions. However the 
behavior is completely different for each of the forces that already have very different predictions of 
Aifi{q = 0,T = 0). SKM yields a nearly constant behavior followed by a strong increase of .4.1^0(9) 
for very large q whereas SLyb decreases to a local minimum at q = qc — 500 MeV. The behavior 
resulted by the use of SKM force shows qualitative agreement with the results by Kaiser within Chiral 
Perturbation Theory For the force SLyb the spin symmetry energy coefficient may decrease 

still more for large values of q, at zero temperature, eventually it may become negative at larger 
densities. The instability associated to ^1,0 < is the one towards a ferromagnetic alignment which 
has been found in several works with several Skyrme forces and relativistic models at higher densities 
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[14, 21, 20, 12, 16, 15|. However this transition is absent in several calculations. The most well 



known calculations in which the ferromagnetic alignment is not found are those ones based on NN 



interactions with different methods 11^, Qy, 18, 0|. However there are particular Skyrme forces which 



do not provide this ferromagnetic phase for nuclear matter: those parametrizations with the inclusion 
of NN tensor Skyrme- type force by Liu et al |72| or using SLyb at low momentum as seen in figure 



3, for higher densities and n-p asymmetries - seen in the second of the references |21|. The effect 
of the momentum dependence, however, is the decrease of ^i,o(^)- Whereas the functional density 
formalism with Skyrme forces and the relativistic (mean-field) models with nucleon-mesons couplings 
are effective models for the nuclear many body problem the NN based calculations are subject to 
approximative methods which may not capture all the relevant degrees of freedom appropriatedly in 
each part of the nuclear phase diagram. At finite temperatures Aifi{q,T) does not vary significantly. 

In the Figure 4 the scalar polarizability, ^o,o as defined in expression ( |3^ ) is plotted. It shows a 
continuous increase with momentum without the turning point at (7 ~ SOOMeV. This parameter, a 
dipolar incompressibility, is proportional to the nuclear matter incompressibility, like it is shown in 
the Appendix. 

In all examples shown above, the increase of temperature is more relevant for larger q and the 
increase of the nuclear temperature always yields larger As,t- Usually it is not expected a large 
variation of the static symmetry energy with the temperature from microscopic calculations in finite 
nuclei ||7^, |7^. Modifications of the chemical potential at high temperatures can lead to stronger 
dependences on T. 

To understand the behavior of the generalized polarizabilities the total densities (A^'^, p"^, M"^), 
as defined above, are shown in figure 5 as functions of q for the parameters of force SLyb. They are 
obtained from the zero-frequence of the generalized Lindhard functions and they generate the behavior 
of expressions (^) seen in figures 1 to 4. Whereas Hq and II2 present a smooth behavior towards to 
zero with the increase of q, the momentum density, |j| oc M, has a dramatic changement at q ^ 2kF. 
This prevents the polarizabilities As,t to become negative for the forces investigated in this work, in 
particular the neutron-proton ^0,1 and spin ^1^0 ones. 

4.2 Other considerations 

From the stability analysis of section 2, the results shown in figures 1 to 4, mainly for the force SLyb, 
suggest that nuclear matter is close to undergo phase transitions around q ~ 500 MeV, i.e., when the 
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exchanged momentum q is nearly twice the momentum at the Fermi surface, kp. The q-dependence 
of the Lindhard function yield Mi{q), as well as Ni{q), pi{q) for oj = 0, which prevents nuclear matter 
to undergo phase transitions. 

This analysis was done for zero energy with As,t only as a function of exchanged momentum. The 
frequence dependence of the polarizabilities was analysed associatedly to the exchanged momentum 



for the dipolar collective motions where zero-sound like excitations were found 11]. In nearly 
symmetric nuclear matter they disappear at temperatures of the order of T ~ 7 MeV ^] (and 
higher temperatures for non zero asymmetries). Their disappearance may occur with the liquid-gas 
phase transition [77, |7^. The increase of the giant dipole isovector resonance width stops so that the 
corresponding energy is probably being used for changing the phase of the system 

These results can also be expected to yield consequences for the Supernovae mechanism and proto- 
neutron or " neutron" stars with their dynamical behavior involving energy and momentum dependence 
of As^t- The symmetry energies contribute, among other ways, by means of ^o,0) -^1,0 (due to the 
coupling to neutrinos) and Aq^i, Ai^i for the different neutron-proton densities and the other related 
effects [f79| , |2^ . The neutronization of a proto-neutron star in the quasi-static phase of the supernova 
can be partially suppressed due to the eventual increase of the symmetry energy coefficient although 
the momentum dependence shown in Figure 1 presents the opposite trend of decreasing ^o.i up 
to g ~ Ikp. This second (dynamical) effect seemingly would facilitate the neutronization and it 
should compete with the former. On the other hand the spin symmetry energy is strongly dependent 
on the used Skyrme interaction. Although the (continuous) increase of ^i^o seems to be rather in 
agreement with other works |l^, |l^ new developments are needed and they include the need of new 
parametrization of effective forces focusing on spin-dependent observables from nuclei and nuclear 
matter. 



5 Summary 

In this paper the nuclear matter symmetry energy terms were investigated as generalized polarizabil- 
ities. Stability conditions with respect to neutron-proton fluctuations were derived being complemen- 
tary to others usually investigated A differential equation for the simultaneous dependence 
of the generalized symmetry energy coefficients on the neutron-proton asymmetry and on the total 
nuclear density was proposed with equation ([TtI) . For this, no considerations about the microscopic 
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reasons for the resulting stability density with a given n-p asymmetry were raised. The stability 
density is, in this case, to be a function of the n-p asymmetry as it should be in a general formu- 
lation. Some solutions for this equation were given. This procedure is interesting for finite nuclei 
as well. These results may be of interest for the investigations of the role of the symmetry energy 
on observables in Radioactive Ions which are being prepared and done mainly for the RIA and GSI 
machines. At different densities the isospin dependence of the symmetry energy may be different from 
the usual one. Finally, within the framework of the linear response of non relativistic nuclear mat- 
ter with Skyrme forces, the g-dependence (exchanged momentum between the components of nuclear 
matter, eg. neutrons and protons) of the coefficients As,t was investigated. For low momenta the 
n-p symmetry energy decreases (linearly and quadratically) until qc — 500 MeV in agreement with 
earlier investigations of the symmetry energy potential ||7l|. In this range of momentum transfer q 
phase transition(s) may take place if other conditions are present, such as different (higher /lower) 
nuclear densities. This indication can be seen in the other symmetry energies, the spin-dependent 
ones, which however depend strongly on the particular Skyrme effective interaction. The results in 
the spin-dependent channels show no defined sign of such ferromagnetic phase for the SLyb force at the 
saturation density. However the decrease of ^i^o('i') with increasing transferred momentum may favor 
such phase transition in different conditions of densities and n-p asymmetries. The scalar coefficient 
•^q,o{q)} the dipolar incompressibility, has continously larger values with the increase of exchanged 
momentum. 
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Appendix: Skyrme force parametrization, functions Vi, relation be- 
tween Koo and Aofl{q = 0) 

In this appendix we exhibit the functions Vi (for expression (34) with parametrization of Skyrme 



forces SKM, SLyb and others ||52|] given by: 

vi2 = toil + xoP.)<5(ri - ra) + I (1 + xiP^) [^(ri - r2)k^ + ^'^^(ri - r2] 

t2 (1 + X2Pa) k'.5(ri - r2)k + 1(1 + X3P.)(ai(pi + pz)^ + a2p")<5(ri - ra), 

o 



(A.l) 



where is the spin exchange operator. The parameters for the forces SLyb and SKM are given 



respectively in references |46, 47 1. 

Prom the hnear response calculation for a time dependent Hartree Fock frame, in the lines discussed 
in |1C, 53], we can write the corresponding functions V^^ and V^^ in each channel for the more general 
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calculation in asymmetric nuclear matter: 
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'2 ' /2''" "'^^ "24 
where pn, Pp and /) are the proton, neutron and total densities of asymmetric nuclear matter, a,b,c 
are the asymmetry coefficients defined in section 4. 

For the longwavelength limit of ^0,0 (^ = 0, P, = 0) the following relation is obtained in terms of 



the incompressililiyt modulus |]ll|: 



^TF + 2Viklpo-%p^+\ 
5 4 



(A.3) 
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They have different relevant terms such as the one proportional to the nucleon kinetic energy at the 
Fermi surface, Tp, and a term from the density dependence of the Skyrme forces proportional to t^. 
This can be seen, in general, by remembering that the calculation of Aofi{p,q = 0) was done with 
the quadratic form for the binding energy in the presence of an external perturbation which induces 
density fluctuations of expression (|3|). It is rewritten below: 



P 



+ e'5p. 



(A.4) 



Terms containing ((J/?)", for n 7^ 2, in were neglected and would correspond to the terms which 
yield the usual i^oo- This more general parametrization will be considered in a forthcoming work. 

The general structure of the zero frequence (real) generalized Lindhard functions II2N at zero 
temperature can be written as: 



no(r = 0) 



M*ki 



, M*kl 

n,(T = o) = ^ 



+ - 1 



2, + 2,qH'^ + -{l-?,cq^){l 



)Ln 



q — 2kp 



q + 2kF 



2k'p 



)Ln 



q — 2kF 



q + 2kF 



(A.5) 
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where ai,bi,Ci,di, Cj dependend on kp and on M* 



Figure captions 

Figure 1 Neutron-proton symmetry energy coefficient ^0,1 = /o/(2n^^) of symmetric nuclear 
matter as a function of the momentum transfer between neutrons and protons, g(MeV), for interactions 
SLyb for T = 0, 4, 7 MeV and SKM (T=0). 

Figure 2 Spin symmetry energy coefficient Ai^ = p/{2Il]f) of symmetric nuclear matter as a 
function of the momentum transfer between neutrons and protons, q (MeV), for interaction SLyb for 
T = 0, 4, 7 MeV and SKM (T=0). 

Figure 3 Spin-isospin symmetry energy coefficient Ai^i = p/{2Il]^^) of symmetric nuclear matter 
as a function of the momentum transfer between neutrons and protons, q (MeV), for interaction SLyb 
for T = 0, 4, 7 MeV and SKM (T=0). 
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Figure 4 Scalar symmetry energy coefficient Aq^q = p/(2n^ ) of symmetric nuclear matter as a 
function of the momentum transfer between neutrons and protons, q (MeV), for interaction SLyb for 
r = 0, 4, 7 MeV and SKM (T=0). 

Figure 5 The densities N, p, M as functions of the transfered momentum between neutrons and 
protons for the force SLyb. They are nearly independent of the force unless for the values of m* and 
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Figure 3 
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Figure 5 
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